The missing asymptotic sector of rotating black-hole spectroscopy 



Uri KeshelEl and Arnon Ben-Meir 
Physics Department, Ben-Gurion University of the Negev, Be'er-Sheva 84105, Israel 

(Dated: July 11, 2012) 

The rotation of a black hole splits its spectrum in two, yet only one sector is known in the 
highly-damped regime. We find the second, at least partly oblate sector, with quasinormal modes 
approaching the total reflection frequencies uj{n S> 1) — mfl — 2niT{n — s), where Q and T are 
the horizon angular velocity and temperature, s is the field spin, and m is an azimuthal eigenvalue. 
Some physical implications are discussed. 
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The natural vibrational modes of a black hole, known 
as quasinormal modes (QNMs) due to their typically fast 
damping, play an important role in the classical study of 
black holes, in the search for a theory of quantum gravity, 
and in the gauge/gravity duality; see [H 0] for reviews. 
The highly-damped modes, in particular, tend to show an 
interesting structure with simple or no dependence upon 
the vibration parameters. It was suggested that such 
an asymptotic QNM corresponds directly to a quantum 
transition, and for a Schwarzschild black hole may be in- 
terpreted as an elementary change in the black hole area 
Q . The subsequent flurry of black hole spectroscopy re- 
search produced many results, but for the most part, the 
asymptotic QNMs defied a simple interpretation along 
the lines of the correspondence principle [l| . 

The asymptotic QNMs of rotating black holes were 
computed numerically 0] and analytically 0. These, 
and the related total reflection and total transmission 
modes, show a curious structure, which provides inter- 
esting hints at the microscopic description of rotating 
black holes in the asymptotic regime Q . Curiously, while 
the rotation is known to split the spectrum into two dis- 
tinct branches of QNMs at low frequencies only one 
asymptotic branch has been reported thus far, and was 
identified with equatorial vibrations @. It is natural to 
ask what happens at asymptotic frequencies to the sec- 
ond QNM branch, and to polar excitations in general. 
The question is further motivated by some contradictory 
results in the literature, and because the very existence 
of a second asymptotic sector was largely overlooked. 

We compute the QNMs of Kerr black holes up to high 
overtones, in order to uncover the fate of the elusive sec- 
ond branch. The transmission-reflection problem is then 
analyzed in the oblate limit, relevant to the new sector. 

Teukolsky's equations. — Consider an electrically neu- 
tral black hole of mass M and angular momentum J = 
aM. Linear, massless field perturbations of the black hole 
are described by Teukolsky's equation [8]. Decomposing 
the perturbation's wave-function into radial and angular 
parts, sMx) = e'('"'*-'"*),5/™(/i = cos 0),i?i,„(r), leads 
to separate angular and radial equations linked by a cou- 
pling constant s^im, which is closely related to the Carter 
constant [§], as we argue below. Here, x = (i, r, 0, 0) are 



Boyer-Lindquist coordinates, and m are angular, az- 
imuthal harmonic indices. The spin- weight parameter s 
gives the spin of the field, specializing the analysis to 
gravitational (s = —2), electromagnetic (s — —1), scalar 
(s = 0), or two-component fermion (s = —1/2) fields. 
We use Planck units where G = c = fcs = fi.= l, and 
shall henceforth omit the indices s^j_, m when possible. 
Teukolsky's angular equation is 
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The regularity of S at the poles /i = ±1 picks out a dis- 
crete set of solutions, known as spin- weighted spheroidal 
wave functions [lol and references therein] and corre- 
sponding eigenvalues A{uj). For s — 0, these functions 
reduce to the familiar spheroidal wave functions 
The radial equation can be written in the form [5| 







(2) 



where A 



-2Mr + a'^ vanishes at both horizons, r± 



M ± (M2 - a2)i/2. Here we defined (coV)^ 
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and q = 



(3) 

E 2i[rA-Af(r2-a2)]a;-(Af-r)[2iaTO+s(M-r)], 
A + s-m'^. The QNMs are solutions of Eqs. dTj) 



and ([2]) for physical boundary conditions of purely outgo- 
ing waves at both spatial infinity and the event horizon, 
i.e. crossing into the black hole. Equivalent ly, these so- 
lutions correspond to transmission and reflection by the 
black hole when the incident wave is negligible. 

For given black hole (M, J) and field {s,l,m) param- 
eters, the above constraints pick an infinite, discrete set 
of QNM solutions. These modes are labelled by an in- 
dex n G N, and are specified by their complex frequency 
Lu — ujn + iuji and coupling constant A, where w/ < 
(time decay) diverges as n — >■ cx). They show a complex- 
conjugate pairing symmetry, where for each {n,/,m}- 
QNM there is an {n,l, — mj-mode that satisfies [7] 



n,l, — ra 



and ^n,/,m — A 
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We thus consider only m > 0, so each mode is unique. 

In addition to the pairing symmetry, the rotation of 
the black hole introduces a Zeeman-like splitting of the 
QNMs, such that for given parameters {n,l,m > 0} 
there are actually two modes, each belonging to a dis- 
tinct QNM branch Q, which we label + (new branch) 
and — (known branch). Each such QNM pair coalesces 
into a single Schwarzschild mode as a — t- 0. 

Series solution. — For purely imaginary large cj, S is 
thought to be of the prolate-type, such that pj| 



A'-P'-^ = tAiauj + Ao + 0{\auj\-^) 



(5) 



where Ai = 2L + 1, L = I ~ max(|m|, |s|) and 
is a constant; Aq = - L{L -I- l)/2 - 3/4 for s = 0. It 
was previously argued that this applies to all QNMs in 
the highly-damped regime [3]. However, this assumption 
was found [1, Q to yield only one, equatorial [g*] branch 
of QNMs, which therefore questions the validity of the 
assertion. To search for all QNMs, here we make no as- 
sumption regarding A, but rather solve both the angular 
equation ([1]) and the radial equation ^ simultaneously. 

QNM frequencies are usually computed by series ex- 
pansions of R and S, each leading to its own three- 
term recurrence relation. These commonly-used recur- 
sions were derived in and will not be reproduced here; 
a simplified form of the radial recursion may be found in 
[l^ . The QNMs are determined by requiring that both 
recursions simultaneously converge, leading to an infinite 
set of QNM parameter pairs {a;„. An}- 

In the highly-damped regime, n ^ 1, the series com- 
putation becomes increasingly challenging numerically. 
In order to find such QNMs, including both branches, 
we start with pairs of easily computed ^13i] Schwarzschild 
modes related by a;+ = —uj~, defined such that 3fJ(a;+) > 

(new branch) and 5i(a;~) < 0. Then we gradually 
increase a from to its extremal value, in sufficiently 
small steps to allow our root-finding algorithm to con- 
verge. The results are illustrated in Figures [T] and [21 

Figure [1] depicts A for gravitational modes with —s~ 

1 = m = 2 (henceforth: fundamental modes). In the 
Schwarzschild limit, A{a = 0) = l{l + I) - s{s + 1) for 
both branches Q ; A~^ , A~ remain comparable at small 
values of the rotation parameter a. However, beyond a 
certain threshold a^., the two branches show markedly 
different behaviors. While the known branch asymptotes 
to the familiar prolate scaling of Eq. ([5]), the new branch 
asymptotes to the oblate limit, given by [3] 

= ~{auj)^ + 2qaLU-iq^-m.^ + 2s + l) + Oi\auj\-'^) , 

(6) 

where q was derived in flO], and in our case reduces to 

fl - TO -I- 2(; + s) iim>l + 2s; 

q= < \ ' ^ (J) 

I 1 + Z — mod(^ + TO, 2) otherwise. 

The split point becomes gradually smaller at higher n; 



comparing Eqs. ([5]) and ([6]) suggests that Og ~ l^i/w/j, 
such that fls -H- in the highly-damped limit. 
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FIG. 1: Absolute value of the coupling constant A for funda- 
mental modes in the missing branch (circles) and in the known 
branch (diamonds), for n = 30 (small symbols), 40 (medium), 
and 50 (large). The corresponding oblate (solid curves; Eq. [6]) 
and prolate (dashed curves; Eq. [S] with Ao = 9) formulae fit 
well above an increasingly smaller threshold as{n). 

Prolate and oblate behaviors are known to interchange 
in the complex uj plane, along branch cuts emanating 
from points where different eigenvalues A coalesce [e.g., 
T^, and references therein] . Our results suggest that such 
branch cuts are present near the imaginary axis for large 
1^1, at least for the fundamental mode (but see [lo|). 
Note that both oblate and prolate asymptotic expansions 



are highly degenerate, A 



(ob) 



A'i'^l when < / - 



2s 



is even, and (at least when s — 0) Am^^ ~ A\ _ 

Figure [5] shows the real (figure body) and imaginary 
(inset) parts of the QNM frequencies of various modes. 
Consider first the new branch of fundamental modes with 
TO > 0, depicted as empty circles of size corresponding to 
n. It quickly converges to ujr = mO. (solid curve), where 
Q, — Ana/ A and A = 47r(r^ -f a^) are the angular velocity 
and surface area of the (outer) horizon. Non- fundamental 
TO > modes do not converge as rapidly with n, but 
are seen to approach or oscillate around this value. All 
modes with to = (filled symbols) rapidly converge to 
ujR = 0. The inset depicts -lui/2t:T - [-uji/2ttT\ for 
various modes, where T = (r+ — r_ ) is the Bckcnstein- 
Hawking temperature. It suggests that for s G {—2, — 1}, 
LOi — !> —2ttT times an integer; for s G {0, —1/2} we can- 
not show such convergence. 

All new branch to > modes substantially differ from 
the known QNM branch, for which ujn is depicted as a 
dashed curve and w/ is nearly a multiple of 2TrT{a = 0) 
!5] . The results suggest that the QNM frequencies of the 
new branch tend to the limit (solid curve in Figure [5]) 



uj{n oo) — mfl — 2'KiT{n — s) . 



(8) 



In order to derive this result, we compute the QNMs un- 
der the assumption that A is given by the oblate Eq. ^ . 
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FIG. 2: The real frequency ljr of missing branch modes with 
s — ~2 (circles), s — —1 (squares), s — (right triangles), 
and s — —1/2 (up triangles), for both m > (empty symbols) 
and m — (filled symbols). Large symbols correspond to 
n = 60. For s — —2, we also show results for n = 40 (medium 
symbols) and n = 20 (small symbols), in order to illustrate the 
fast convergence of new branch modes onto ujr = mil (solid 
curve) and the slower convergence of known branch modes 
(diamonds) onto the known solution (dashed curve). Inset: 
the fractional modulo of the imaginary frequency \ioi\ with 
respect to 2nT (same symbols). 

This yields Eq. numerically and analytically, as we 
show next. However, while oblate behavior is evident for 
the fundamental modes in Figure (TJ numerical conver- 
gence is insufficient to establish the asymptotic angular 
behavior of other modes, as we discuss below. 

Transmission and reflection. — To establish the new 
QNM branch and study the corresponding asymptotic 
sector, we solve the transmission and reflection problem 
in the oblate, highly-damped limit. Our analysis uti- 
lizes the Stokes phenomenon and monodromy matching, 
in resemblance of the prolate study in 1^. However, 
the oblate geometry considered here leads to a different 
complexified geodesic structure, illustrated in Figure [31 

We define a tortoise coordinate z = V dr, such that 

z r as r ^ oo ; z— > — cxd as r r+ , (9) 
and Eq. ([2]) becomes 

where f{z) = V^A^R and V = 3V'{r)^/{4:V^) - 
V"{r)/{2V^). For an incident wave of amplitude I, 

^ |Xe-»"^ + 7e(w)e'"^ as r ^ oo , z -J> oo ; ^^^^ 
I T(w)e~''^^ as r ^> r+ , z — s> — oo , 

where 7" and TZ are called the transmission and reflection 
amplitudes when 1=1. Here, we imposed purely outgo- 
ing boundary conditions at r^^ so f is an eigenvector of 



the monodromy matrix with eigenvalue e^'^"'^+ , where 
ei^-^i = e^^-'«/2 . (12) 

For a discussion of the subtleties involved, see @. 

We derive cj, T and TZ in the WKB approximation, by 
evolving / in the complex r-plane along anti-Stokes lines 
defined by 3?(iwz) — 0. These lines, identified below 
as complexified geodesies, are shown as solid curves in 
Figure [3] and emanate from the turning points in which 
V = 0, depicted as circles. The contour C we choose 
is analytically continued from r — cx) to the geodesic 
li, continues from there to I2 after avoiding the turning 
point t_ by a 2tt/3 counterclockwise rotation, to ^3 after 
a clockwise rotation around r+ , to li after a counterclock- 
wise rotation around t^, and finally back to r ^ cx). 
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FIG. 3: Illustration of anti-Stokes (solid curves) and Stokes 
(dashed curves) lines in the complex r = r_R + iri plane. Also 
shown are the horizons (disks) and the turning points (cir- 
cles). Here, a = M/2 and Muji = -500. 

Monodromy matching and the mixing of WKB modes 
imply, in resemblance of the prolate analysis [6] , that 

r{uj) = + ^7^(w); r(-w) = X(-a;) = ^7^(-w); (13) 
and e-4'^"'^7^(a;) = 7^(a;) + [iX(a;) - 7^(w)]e-2'"'^ , (14) 

where the integral in 5 = z(t+)_-z(t_) is evaluated along 
C in a single Riemann sheet [25|. In the highly-damped 
limit, e"^*'^'' 0, so Eq. yields the frequency condi- 
tion Eq. ([5]), as anticipated above. Equation ([T3)) guar- 
antees flux conservation, T(w)T(— w) + TZ{uj)TZ{—uj) = 
I{uj)I{—u!), but does not entirely fix T{io) and TI{lu). In 
particular, these equations admit QNMs with I{uj) — 
and total reflection modes with T{uj) — 0, with the same 
asymptotic frequencies. For total reflection, Eq. ([5]) is 
precise even for finite w, as in the prolate case [6|. 

Discussion. — After pointing out that previous stud- 
ies of rotating black hole spectroscopy in the highly- 
damped regime were limited to only one out of two 
sectors, we provide numerical evidence that the miss- 
ing QNMs asymptotically approach the total-reflection 
branch Eq. ([5]). This is shown in Figure for all spins. 
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and is particularly striking for the converged fundamental 
mode. Our results indicate that the new sector involves 
oblate angular wave functions with eigenvalues A given 
by Eq. ([5]), at least for the fundamental mode, implying 
a branch cut near the imaginary w axis. Such behavior is 
seen in Figure [TJ and is consistent with our numeric and 
analytic computations (Eqs. ITMMI) . which show that 
oblate behavior directly yields Eq. ([8]). 

However, the prolate analysis [5| , which found only one 
QNM branch in the highly damped limit, breaks down for 
the frequencies of the new branch, where the two mon- 
odromy eigenvalues in Eq. (IT2]) become degenerate Q. 
This is not the case for our oblate analysis, where the 
QNMs gradually approach Eq. ^ as e"^*"'' 0. There- 
fore, we cannot rule out the possibility that for some 
non-fundamental modes, the new branch is prolate. Our 
numeric results cannot reliably resolve this issue, due to 
the difficulty of computing high overtones. 

Consider the angular behavior of the new branch in 
the highly damped, oblate case. Equation ([6]) indicates 
that in the Teukolsky Eq. ([T]), the real part of the square 
brackets on the right hand side is large and negative 
away from the poles, so arg(S') oscillates wildly there. 
This suggests that S is exponentially small away from the 
poles, which is indeed the known behavior for s = [e.g., 
Il6j . and appears to hold numerically for all s. At high 
overtones, the oblate branch thus involves polar modes, 
in contrast to the prolate branch associated with equato- 
rial modes; these roles are reversed at large real uj. Note 
that it is the new branch Eq. ([8|) which at small n ap- 
proaches superradiance in the extremal limit. 

Consider the geodesic motion of a massless particle 
with angular momentum = m, complex energy E — uj, 
and Carter's (fourth) constant of motion ^ fixed at some 
value qc. The covariant radial momentum is then given 
[13] precisely by pr in Eq. jS]), provided that qc = q. 
Based on both prolate and oblate results, we thus pro- 
pose that Carter's constant is to be identified with A, up 
to an energy-independent term. The potential ujV can 
be understood as a generalized momentum p which is the 
quadratic mean of pr and a spin term sps. This justifies 
the geodesic interpretation of Figure (Sj for a more de- 
tailed discussion, see [6]. As in the prolate case, the QNM 
condition ^ can now be written as a Bohr-Sommerfeld 
equation, 2 §^pdr = hn, where h is Planck's constant. 

The real part of the new QNM branch agrees with the 
first law of black hole thermodynamics, AM = TAS + 
flAJ, if the black hole undergoes an angular momentum 
change A J — m and an entropy change AS = |18|]. In 
this sense, the new QNMs admit a simple interpretation 
as reversible angular momentum transitions of the black 
hole. Modes with m — entail no change in M or J; such 
lor = modes could correspond to forbidden quantum 
transitions Interestingly, while attempts to infer 

the area spectrum from the QNM level spacing [e.g., 



yield Bekenstein's 2lj result A^ = Stt only in the slow 
rotation limit using the known QNM branch [23 | , the new 
branch gives Stt for arbitrary a. 

In hindsight, several clues to the nature of the miss- 
ing branch may be found in the literature. The behavior 
uju = 2il was noticed |23|] for the fundamental mode, 
and associated with reversible black hole transitions [3| , 
but this was considered an intermediate or unusual case, 
and was later contradicted [e.g., 14. \ lp\ . A result similar 
to Eq. dH) was postulated in Ref. [2j| based on a corre- 
sponding ~ [oj[~^ drop in the recursion coefficients in the 
oblate limit, which may in principle truncate the expan- 
sion series. Finally, an expression reminiscent of Eq. 
was found at intermediate asymptotic frequencies in [121 ■ 

We thank S. Hod for encouragement and helpful dis- 
cussions and suggestions. 
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